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PsJ ' Abstract 

C ' Starting from the canonical phase space for discretised (4d) BF-theory, we implement a canon- 

ical version of the simplicity constraints and construct phase spaces for simplicial geometries. Our 
construction allows us to study the connection between different versions of Regge calculus and 
approaches using connection variables, such as loop quantum gravity. We find that on a fixed 
triangulation the (gauge invariant) phase space associated to loop quantum gravity is genuinely 
'~~i. larger than the one for length and even area Regge calculus. Rather, it corresponds to the phase 

^ji' space of area-angle Regge calculus, as defined in [1] (prior to the imposition of gluing constraints, 

I , which ensure the metricity of the triangulation). Finally, we show that for a subclass of triangula- 

^jq' tions one can construct first class Hamiltonian and Diffeomorphism constraints leading to fiat 4d 

I— I. space-times. 

^ ■ 1 Introduction 

o: 

OO ■ In many approaches to quantum gravity, such as loop quantum gravity [2j, quantum Regge calculus 

^^ . [U Hj and causal dynamical triangulations [5], discrete models of space-time appear in one form 

t^ I or another. One advantage of these discrete models is that they have the potential to provide a 

lyQ • regularisation of uv divergences [BIE]. On the other hand, discretisation very much complicates the 

^^ . role of diffeomorphisms, the underlying symmetry of continuous general relativity. Moreover, as we 

^ I will show in the course of the paper, even on the kinematical level, discretisation can lead to the 

appearance of non-metric configurations which are absent in the classical and continuum phase space 
description. 
5^ I These complications arise in both the path integral quantisation and canonical quantisation. Spin 

foam models and quantum Regge calculus involve a path integral over discretised geometries. Were 
diffeomorphisms to act on the configurations within these path integrals, divergences might appear, 
caused by an integration over the diffeomorphism gauge orbits [8j. Thus, it is important to under- 
stand how diffeomorphisms act on discrete manifolds. In a canonical quantisation based on discretised 
structures, the typical problem is to obtain a consistent time evolution. In the continuum theory, 
time evolution is generated by the Hamiltonian and diffeomorphism constraints, which additionally 
generate 4d diffeomorphisms (on-shell). The group property of these diffeomorphisms is reflected in 
the fact that the Hamiltonian and diffeomorphism constraints are first class. In other words, in a phase 
space description, the Poisson brackets of these constraints vanish on the constraint hypersurface. Un- 
fortunately, this first class property is typically lost [9l|Tni[II] if the spatial geometries are discretised. 
One could argue that discretisation breaks diffeomorphism symmetry [12], and so we do not need to 
care about diffeomorphisms at the discretised level. Nevertheless, diffeomorphism symmetry has to be 



X 



restored in the continuum limit, as it is a symmetry of the continuum theory. Therefore, one would 
expect at least some approximate notion of diffeomorphism symmetry in the discretised theory. More- 
over, trying to preserve some kind of diffeomorphism symmetry (that is, trying to preserve the Dirac 
algebra of Hamiltonian and diffeomorphism constraints) could help to resolve quantisation ambiguities 
[T3j in the Hamiltonian constraints. 

Although there is a considerable amount of work addressing these questions in different frameworks 
[TT| [Ml [151 HSl llTl [HI [iHl [20], in our view the issue is not yet settled, even at the classical level. 
One result in this paper is that one can obtain a phase space description for (3 + 1) -dimensional 
discrete geometries with the following properties. One can impose a consistent dynamics, at least 
on a certain class of special triangulations, which leads to flat 4d geometries and upon which there 
is a well-defined action of the diffeomorphism constraints. Since flat geometries are a subsector of 
the space of all solutions (of the Regge equations, for instance), we expect at least some remnant of 
the diffeomorphism symmetry from the full description of the theory. In fact, given a simplicial 4d 
manifold that satisfies the Regge equations of motion, one can always subdivide a 4-simplex of this 
triangulation by placing a new vertex into this 4-simplex such that the geometry of the 4-simplex is 
still flat. In this way, one obtains four gauge degrees of freedom, as the geometry of the triangulation 
is not changed and the vertex can be freely placed inside the 4-simplex. 

Our results are similar to the work [TTj by Waelbroeck and Zapata which defines 'topological 
gravity' as a subsector of BF-theory. This subsector has 'flat dynamics' and they impose conditions 
which allow for a geometric interpretation of the bivector field E and the curvature field F. These 
conditions are a canonical version of the simplicity constraints appearing in the Plebanski formulation 
|21| of general relativity. We will differ, however, in the methods and in certain parts of the conclusions. 
In particular, we will perform a full reduction of the BF-theory phase space to a phase space describing 
geometric configurations. This phase space describes 'kinematical configurations', and can be taken 
as a starting point to define a dynamics leading to the full set of (Regge) gravity solutions. 

Another motivation of this work is to explore connections among different quantum gravity ap- 
proaches, mainly among spin foam models, canonical loop quantum gravity (LQG) and the different 
versions of Regge calculus, namely length Regge calculus |22j, area Regge calculus [23J and area-angle 
Regge calculus [Ij . See for instance [241 [25 ] for attempts to connect the Ashtekar variables on which 
LQG is based, with Regge calculus. We will find a correspondence between these theories only at the 
kinematical (and classical) level as we do not discuss the full implementation of the dynamics here. 
Nevertheless, if there is a definition of the diffeomorphism and Hamiltonian constraints available in 
one of these versions the methods presented here allow one to translate these to the other versions. 

The definition of spin foam models [26] usually starts with the Plebanski action, which can be seen 
as a BF-theory action plus additional terms involving the simplicity constraints. These constraints 
impose that the bivector field appearing in this action can be written as (the Hodge dual of) a wedge 
product of tetrad fields. A considerable amount of current work on spin foam models [271 [Ml [SH [30] 
discusses the question of how to properly impose the simplicity constraints. This question is crucial 
as BF-theory in itself is a topological theory, that is, it does not have local degrees of freedom. 
Only if the variations of the bivector field are restricted by the simplicity constraints does one obtain 
local degrees of freedom. Similarly in area(-angle) Regge calculus one uses areas (and angles) as 
fundamental variables. These also need to be constrained to come from a consistent assignment of 
length variables to the edges of the triangulation, in order to obtain a theory with local degrees of 
freedom. 

We therefore follow [IT] and start with the phase space of discrete SO {4) BF-theory. In addition 
to the canonical structure defined in |17j , we introduce the Barbero-Immirzi parameter following |31j , 



to keep in line with the starting point for many of the spin foam models (and to a lesser extent loop 
quantum gravity). On this phase space we will impose the so-called simplicity constraints ensuring 
the geometricity of the configurations, i.e. that the bivector field follows from the dual of the wedge 
product of vectors associated to the edges of the triangulation. The aim is to obtain a reduced phase 
space with respect to these constraints and hence a phase space describing simplicial 3d geometries 
embedded into 4d geometries. To facilitate the geometric interpretation we will work with 5'0(4)- 
gauge invariant quantities, that is areas, 3d and 4d dihedral angles and length variables. This also 
allows the connection to the different forms of Regge calculus. 

As we will see, the simplicity constraints can be subdivided into three sets which can be separately 
dealt with. We will obtain the phase space corresponding to loop quantum gravity restricted to a 
fixed triangulatioro already after we have imposed the first set of these constraints. Thus, in this 
phase space the simplicity constraints are not fully implemented and therefore degrees of freedom are 
included which do not correspond to (lattice) metric degrees of freedom. The question arises of how 
spin foam models with proper implementation of the simplicity constraints can be connected with 
canonical loop quantum gravity, where the canonical version of these constraints are not implemented 
fully. 

Of the three subsets of simplicity constraints mentioned before, the first two sets of simplicity 
constraints are second class in themselves, however the last set - giving constraints between the 
areas only - is not. Here one either reduces by these first class constraints or introduces further 
constraints conjugated to these area constraints. These further constraints might arise as secondary 
constraints from the primary simplicity constraints, since in the continuum |32j there are secondary 
constraints coming from the Poisson brackets of the simplicity constraints with the Hamiltonian and 
diffeomorphism constraints. 

Such left-over simplicity constraints do not appear in (2+l)-dimensional gravity. Thus, also the 
phase space associated to loop quantum gravity restricted to a fixed triangulation corresponds exactly 
to the one from length Regge calculus [33]. The reason why the issue of defining a phase space 
for (3 + 1) simplicial triangulations is more involved than in (2 + 1) dimensions is the following: 
the natural configuration variables for a simplicial manifold are the length variables associated to 
the edges of the triangulation. These length variables specify completely the internal geometry of 
the triangulation and are moreover unconstrained (apart from the triangle inequalities and higher 
dimensional generalisations which we will ignore). The natural variables to encode the extrinsic 
curvature are the 4d dihedral angles between neighbouring tetrahedra and therefore based on the 
triangles of the triangulation. Now in general for non-degenerate simplicial 3d (and 4d) manifolds 
the number of triangles is greater than the number of edges. The problem of left-over simplicity 
constraints is therefore rooted in the discretisation by triangulation, which leads to a mismatch in the 
number of variables, which are, however, in one-to-one correspondence in the continuum. 

The structure of the paper is as follows. In Section [21 we shall define the canonical variables and 
Poisson structures for a discrete BF-theory. In the subsequent section, we shall develop a complete set 
of simplicity constraints for the discrete phase space. We prove that these constraints are sufficient to 
ensure the geometricity of the 3d triangulation in Section [H We introduce gauge-invariant variables 
in Section [SJ which will provide the necessary parameters for the Gaufi reduced phase space. We 
consider the explicit case of the boundary of a 4-simplex, and through Section [6] we enforce the 



^With this phase space we mean the following: we restrict the Hilbert space of loop quantum gravity to a fixed 
triangulation and consider the corresponding classical phase space associated to this triangulation. This phase space 
is different from the continuum phase space, i.e. the Ashtekar connection and conjugated densitised triad fields, upon 
which the LQG quantisation is based. On this continuum phase space, simplicity constraints do not appear. 



GauB and simplicity constraints to arrive at the reduced phase space of geometric configurations. In 
particular, we divide the examination of the simplicity constraints into two subsections. The first, 
subsection 16. 11 deals with the constraints implementing the equality of left- and right-handed sectors. 
The second, subsection 16.21 reduces by the constraints ensuring that the edge lengths are independent 
of the tetrahedron in which they are calculated. In Section [Tj we expand our analysis to more general 
configurations, while in the penultimate section, we consider implementing the dynamical constraints 
on the phase space. Finally, we recapitulate and conclude in Section [9l 

2 Canonical structures for 4D BF— theory 

In this section we will collect some basic material and notation on the canonical formulation of discre- 
tised BF theory as this will provide us with the phase space on which we will implement the simplicity 
constraints. Additional material can be found in Appendices Rl and iBl 

Canonical structures for continuum BF-theory are dealt with extensively in [31] . We shall develop 
analogous structures for the discrete setting. The action for 4d BF-theory with an extra 'topological' 
term is 

Sm[E,A]= f (e + -*e] AF= f Ea(f+-*f], (2.1) 
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where A^ is a smooth manifold, E is an so(4)-valued 2-form and A is the so(4)-valued 1-form. They 
are referred to as the bivector and connection fields, respectively. The curvature of this connection is 
F = dA + A f\ A. We shall refer to the coupling constant 7, as the Barbero-Immirzi parameter even 
in the non-gravitational theory. 

It is straightforward to see that the connection A is canonically conjugate to the field 11 := E-\--*E. 
In fact, we shall wish to deal mainly with the bivector field which we may express in terms of the 
canonical momentum as 

7^ f^ ^ 



^=yeT(^n--*nj. (2.2) 

As is well known, the 50 (4) Lie algebra splits into selfdual and anti-selfdual sectors. We may make this 
distinction exphcit by introducing selfdual and anti-selfdual fields E± = ^{E±*E), A± = ^{A±*A)]^ 
which results in 

S[E±,A±] = (1 + -) f E+AF++(l--j /^_AF_. (2.3) 

Note, the special case 7 = ±1. In this scenario, the action contains only the selfdual fields or the 
anti-selfdual fields, respectively. Thus, only half the number of variables are present in comparison 
to the case of general 7. In this case, the (anti-)selfdual connection A± is directly conjugate to the 
(anti-) selfdual field E±. 

Now, we shall repeat this process in the discrete setting. We describe this undertaking in Appendix 
IBI but we shall reiterate some of the pertinent details here. We pass to the Hamiltonian formalism 
and replace the continuum 3d hypersurface by a 3d simplicial complex. Likewise, we define fields on 
sub-simp lices of this discrete manifold. The discrete analogue of the action is 



S[Uij,Mij,N„Ni]= [ dt 



{ij} e i 



it)- (2.4) 



^This implies F± = F[A±] 



where we label the tetrahedra by the index i, so that one can pick out an oriented triangle {ij} by 
the pair of tetrahedra that share it. Moreover, e denotes an edge of the siniplicial complex. The 
canonical variables are associated to the triangles {ij} where (Ilij,Mij) are the discrete counterparts 
of the canonical variables (11,^). Mij is a rotation between the frames of reference in the adjacent 
tetrahedra, and IIjj = Eij + -*Eij. The latter two terms in (12. 4p are the flatness and Gaufi constraints. 
We discuss their significance later on, but they are enforced by Lagrange multipliers N^, and Ni, 
respectively. The Poisson brackets between the bivectors Eij and the holonomies M^j are then given 
by 

{Ef,, M,r } = :^ {5^^' - '-e^^') {S,k5,, C^'^^M,r + Su6,, C^'^^M.^ ) , (2.5) 



{Ei,E^,} = 5,,6,,^L_c^BC ( ^cc _ l^cc ) ec; . (2.6) 



t 
'^^ — 1 \^ ^ 

Between the canonical variables we have two kinds of relations which are consistent with the 
Poisson brackets (|2.5p . Firstly, Mji is the rotation inverse to Mij and secondly, we have the relation 
Ej- = —M^Ef- which relates the bivector as seen in both tetrahedra adjacent to the triangle {ij}. 

In this setting, we may also perform the decomposition into selfdual and anti-selfdual parts. The 
projected fields and structure constants are E^^^-, M^^ and C^^'-' respectively, which we detail in 
Appendix O In summary, we obtain the following symplectic structures (assuming i = k and j = I) 

{EtMi''} = ^-C^^^Mf^, 

7 lb 1 

{ElEi} = ^Ci^<'Eg. (2.7) 

7 ± i 

The constraints in terms of {Eij, Mij) are the discrete flatness and Gaufi constraints 

pA ^ qAbc (^Mi.Mjk . . . M^if^ , (2.8) 

Gf = Y.^i. (2.9) 

k 

where the product on the right hand side of ()2.8p is the holonomy around the edge e G A, and the sum 
on the right hand side of (|2.9p imposes closure of the four bivectors associated to the tetrahedron {i}. 
Both Gaufi constraints and flatness constraints form a first class constraint set. The Gaufi constraints 
generate 50(4) gauge transformations 

Ef^ ^ Af^'E^, (2.10) 

for 50(4) matrices Aj,Aj associated to the tetrahedra {i} and {j} respectively. The flatness con- 
straints generate the translation symmetry 

M,f ^ M^^, 

E^ ^ Ef,+ Y. ^ei, (2.11) 

em{ij} 

where the Agj are gauge parameters associated to the edges of the triangulation. They are expressed 
in the reference system of the tetrahedron {i} and we sum over the three edges of the triangle {ij}, 
see Appendix |F] for more details on the action of the flatness constraints. 



3 Canonical simplicity constraints 

In this section, we will introduce a canonical version of the simplicity constraints, which will impose 
the metricity of the bivector fields. That is, for phase space points satisfying these constraints, the 
bivectors can be written as (the Hodge duals of) wedge products of vectors that can be consistently 
associated to the edges of the triangulation. 

The action (|2.ip describes a topological field theory, namely BF-theoiy. Adding to this action 
another term, known as the simplicity term 

S[E,(t>Ump = -l [(I>^''E'^AE'' , (3.1) 



one arrives at an action which describes gravity (and another topological theory) j2Hl26j. 4> = (j)^ 
is the Lagrange multiplier which ensures that the bivector field arises as E = ib*e A e or £' = ite A e. 
The first pair of solutions are known to give rise to the gravitational sector of the theory, while the 
second pair give rise to the topological sector. 

As one might expect, we wish to have a comparable set of simplicity conditions for our discrete 
theory. But before we proceed, let us note that this action could be discretised using a 4d triangulation. 
But to perform a canonical analysis of the discretised action is in general quite difficult. A canonical 
analysis of the continuum action has only appeared recently jSO]. One of the problems that a discrete 
approach faces is to identify spatial slices in a general 4d triangulation. Generically, the number of 
variables associated to different spatial slices will vary. One could remedy this situation by considering 
a special triangulation, for instance a triangulation that fits into a hypercubical lattice. One could 
then restrict the class of spatial slices to those which correspond to 'planes' in this hypercubical lattice. 
This would, however, introduce a preferred slicing. 

As can be seen in Section [2l we will follow a different route here and discretise only the spatial hy- 
persurface. In fact, we impose constraints on the canonical phase space of discretised BF~theoiy that 
ensure the geometricity of the configuration. The aim will be to reduce the phase space with respect to 
these constraints and the GauB constraint generating 5*0(4) gauge transformations. The phase space 
obtained in this manner should correspond to the one describing true geometric triangulations. As the 
conditions for such a phase space, we will ask that every phase space point allows one to determine 
(at least locally) a set of consistent edge-lengths for the associated triangulation. Furthermore, we 
require a non-degenerate symplectic form on this phase space. It is on this reduced phase space that 
one could analyse the effect of instituting the diffeomorphism and Hamiltonian constraints. 

We will start to discuss such geometricity conditions for a general 3d triangulation. These con- 
straints will be highly redundant, however, and therefore it is difficult to obtain the reduced phase 
space in the general case. As a result, it will turn out to be enlightening to first consider a simple 
closed 3d triangulation, namely the 3d boundary of a 4-simplex. Later on, we will consider more 
complicated triangulations. 

As a first set of constraints, for every tetrahedron we have the discrete equivalent of the GauB con- 
straints, ensuring that the tetrahedra are closed 



E^: 



,j = 0, (3.2) 



where the sum is over the four tetrahedra {j} adjacent to the tetrahedron {i}. Furthermore, we will 
introduce the well known diagonal simplicity and cross simplicity constraints 

Sij := e E^jE^j = , Sijk ■= e E^^E^^ = , (3.3) 



among the bivectors belonging to one tetrahedron {i}. (These constraints have been used also in 
|17|.) The diagonal simplicity and cross simplicity constraints ensure that every bivector E spans 
only a two-dimensional subspace and that every pair of bivectors from one tetrahedron spans only a 
three-dimensional subspace. 

We can reformulate the simplicity and cross simplicity conditions (|3.3p using the splitting of the 
fields into selfdual and anti-selfdual terms. We define 

Aij± = Eij± ■ Eij± , Aijkj^ = Eij± ■ Eik± , (3.4) 

where Ti ■ T2 = T(^T2- Then, we can rewrite the constraints as 

that is the 3-geometry as defined in the selfdual sector coincides with the three-geometry as defined 
in the anti-selfdual sector, see also |34| . 

These conditions are not sufficient. We have also to ensure that all the bivectors meeting at an 
edge span only a three-dimensional space, as every bivector should be normal to the edge in common. 
Therefore, we add the edge simplicity constraintqfl 

a,ki = e''''EiM,fEfi, (3.6) 

where E;^ and E?i are bivectors associated to two triangles meeting at an edge and belonging to 
neighbouring tetrahedra {i} and {j}. Thus, we do not include all pairings of triangles meeting at 
an edge. If more than five tetrahedra share an edge, it might happen, that two triangles meet in 
an edge which are not from neighbouring tetrahedra. Nevertheless, (by transitivity) the constraints 
(j3.4p are sufficient. Note that the constraints ()3.6p involve - in contrast to the covariant simplicity 
constraints ()3.ip - the connection variables. As we will see later on this set of constraints implies the 
gluing conditions (j6.14p . which are not realised in a phase space associated to loop quantum gravity 
restricted to a fixed triangulation. 

4 Simplicity constraints ensure geometricity 

Here we will show that the simplicity, cross simplicity and edge simplicity constraints ensure that the 
bivectors come from (duals of) wedge products of vectors associated to the edges of the triangulation. 

Note that the diagonal simplicity, cross simplicity and edge simplicity constraints are left invariant 
under the transformation E — )■ *E of the bivectors, which exchanges the gravitational sector with the 
topological sector. Hence these constraints do not distinguish between the sectors; see, however, the 
discussion in Appendix [Pl 

Consider two bivectors, Eij and E^., associated to the tetrahedron {i}. The diagonal simplicity 
and cross simplicity constraints ensure [35] that the bivectors are simple and have a common factor, 
that is 

Eij = u Av , Eik = u Aw. (4.1) 



^This set of constraints can be replaced by an alternative set proposed in [17]. As discussed in Appendix iDl this 
alternative set distinguishes between the gravitational and the topological sector, in contrast to (|3.6|l . 



7'' = 


= Cj s At , 


*Eik = 


= CktAu, 


*Eii = 


= QU A s , 



'ij 
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Xij A Ui , 


'ik 
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Xik A rii , 


\l 
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Xii A Ui , 



For three simple bivectors Eij,Eik and En there are two possibihties to ensure that every pair spans 
only a two-dimensional subspace, namely 

or *Eik = Xik A Hi , (4.2) 

*Eii = : 

where s,t, u and Xij,Xik,Xii,ni are vectors while Cj,Ck and q are scalars. 

Note that the two solutions can be mapped to each other by applying the Hodge dual to the 
bivectors. The right hand side of ()4.2p leads to the so-called topological sector, where the dual 
bivectors are equal to the wedge products between the normal to the tetrahedron, ?ij, and the normals 
to the triangles Xij , Xik , xu . The left hand side corresponds to the gravitational sector where the dual 
bivectors are given by wedge products of edge vectors. 

On the gravitational sector, the dual of the fourth bivector Eim associated to the tetrahedron {i} 
satisfying diagonal simplicity and cross simplicity with the other three bivectors can be written in the 
form 

*Ei^ = {u + as)A{f3s + 5t), 

= (3uAs-6tAu + a6sAt. (4.3) 

By rescaling u and t in ()4.2p we can assume w.l.o.g. that in (14. 2p q = 1 and c^ = 1. Multiplying the 
GauB constraint ^^ *£'j„ = = ^^ Ein with riij A riik, n-ik A nn and riu A Uij respectively (where Uij 
is normal to s and t and so on,) we can conclude that 

P = -l, 5 = 1, Cj = -a. (4.4) 

Finally, we can absorb the factors of a into a redefinition of the one-vectors s' := |a|^'^s, t' := 
|ap^t, u' := |a|~^"u, where in the following we will drop the prime. In summary, we obtain for the 
gravitational sector 

*Eij = ±sAt, 

*Eik = tAu, 

*Eu = uAs, ^^-^^ 

*Eim = (u=Fs) A(t-s). 

We shall denote by Vim the vertex of the tetrahedron {i} which is opposite the triangle {im}. If 
(s, t, u) is a triple of positively oriented vectors starting from Vim, then, with the upper choice of sign, 
()4.5p defines the inward pointing bivectors corresponding to the triangles. If we reverse the direction 
of u in the triple {s,t,u), then (j4.5p with the lower choice of sign defines outward pointing bivectors. 
Redefining u' = —u for this case, so that (s, t, u') is again a positive oriented triple of edges starting 
from Vim-, gives *Eij = —s At, *Eik = —t Au' and so on. 

Hence, we obtain the well known result [35] that the diagonal simplicity, cross simplicity and 
GauB constraints allow for four different types of solutions. For the gravitational sector, either the 
dual bivectors or minus the dual bivectors are coming from a tetrahedron. Since applying the Hodge 
dual to the bivectors preserves the constraints, we have also two other types of solution, namely that 
the bivectors or minus the bivectors are coming from a tetrahedron. 

So far we obtained this result for each tetrahedron separately. But we must also ensure that these 
tetrahedra glue to each other consistently. For simplification, we will restrict our considerations to 
the gravitational sector. 



Consider two neighbouring tetrahedra {i} and {j} and three bivectors Eij, Eik, MijEj^i meeting 
at an edge. Because of the previous discussion, we can conclude that 

*Eij = asAt , *MijEji = a's' A t' , 

*Eik = atAu, *MijEjki = a'u' A t' , 

where (s, t, u) is a positively oriented triple, (s', t', u') is negatively oriented and a,a' = i: encode the 
sign ambiguity. 

From the condition (IB.lSp . we have s A t ~ s' At' and hence both s' and t' are linear combinations 
of s and t. The edge simplicity constraint ()3.6p ensures that the four vectors t,t' and u,u' span a 
three-dimensional subspace. If u' is a linear combination of s, t and u all the simplicity relations are 
satisfied. In this case, however, the two tetrahedra span the same three-dimensional subspace, that is 
the dihedral angle between them is vanishing. Therefore let us assume the generic case, that is that 
the inner product between u' and the normal to the tetrahedron {i} is not vanishing. Since t',u',t,u 
span a three-dimensional subspace we can conclude that in this case t' is a linear combination of t,u 
and u' . At the same time, t' is a linear combination of s and t, 

t' = at + Pu + 5u' = a's + j3't . (4.6) 

Taking the inner product with a vector normal to t, u and u' we obtain q' = and hence t' = at. In 
the same way, we can conclude from the edge simplicity at s that s' = j3s and using as At = —a's' A t' 
we have a = —aa'(3~^. So far, we obtained 
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-^), 


*MijEjm' 



—as A t , 

aa't A u' , 

—a^'^au' A s , 

a'{u' + a^^aa's) A {at + a'^aa's) 



(4.7) 



From the (edge) simplicity between *Sjm and *MijEjm' we have 

at + a~^aa's = X{u — s) + /i(t — s) + ^{u' + a~^aa's) . (4.8) 

Again, by taking the inner product with the vector normal to the tetrahedron {i}, we can conclude 
that u = 0. From the inner product with the normal to {j} we have that A = 0. Hence, at + a^^aa's = 
fit — fj,s, from which we can conclude that aa' = —1 and fi = a = a~^ = 1. Hence, in the generic case 
of non-parallel tetrahedra we have proven that the dual bivectors can be consistently written as wedge 
products of edge vectors also for neighbouring tetrahedra. Moreover the choice of orientation is the 
same, i.e. the bivectors are either inward or outward pointing for pairs of neighbouring tetrahedra. 

5 Gauge invariant quantities 

In this section we will construct a complete set of 50(4) gauge invariant quantities. These will simplify 
the implementation of the simplicity constraints and moreover make the connection to the different 
version of Regge calculus obvious. 

Gauge invariance is not broken by projecting down onto the selfdual and anti-selfdual sectors. 
Thus, we shall find it convenient to define gauge invariant quantities for both the left-handed and 



right-handed sectors in parahel. The easiest gauge invariant quantities to construct are the squared 
areas Aij± and the 3d dihedral angles 

, Aijk± Eij± ■ Eik± 

cos (pijk± := = = . (5.1) 

\l Mj±Mk± V -^ij± • -E'iji Eik± • Eik± 

As already noted, the diagonal simplicity and cross simplicity constraints imply that the areas and 
the 3d dihedral angles computed from the left- or right-handed sectors coincide. 

Another class of variables is the 4d dihedral angles. We define bi-vector 'normals' as NA^ := 
{Eij± X Eik±)^ = ±C±.^^Efj^Efl^ and the angle between two of them aQ 

cos eik,ji±:= ^ ^ ^ (5.2) 

\/Nijk± ■ Nijk± Njii± ■ Njii± 

where the result may not only depend on the choice of tetrahedra {i} and {j} but a priori also on the 
choice of triangles {ik} and {jl}. 

As we will see, this is not the situation in the geometric sector. For the case that the three triangles 
in question share the same edge, we can write 

*Et^' = 2ef ef 1 , *Mi,E'^f = -le'^ef^ , (5.3) 

where (ei, 62, 63) is a positively (negatively) oriented triple of vectors whereas (ei, 62, Cg) is a negatively 
(positively) oriented triple. For the normal bivectors N, we obtain 

iV-; = 2el"nf ] , N^ti = 2eS"nf , (5.4) 

where nf = e ele'^e'^ is the positively (negatively) oriented normal to the tetrahedron {i} and 
n°: = e"^^'^e5e2e'3 is the negatively (positively) oriented normal to the tetrahedron (j). Thus, cosOij^ki 
gives the (inner) dihedral angle between the two normals to the two tetrahedra in the case that 
{ij},{ik} and {jl} have an edge in common. Also note that in this case, 9ik,ji does not depend on 
the choice of the triangles Eik, Eji as long as they share the same edge with Eij. 

To simplify the expression (j5.2p for the 4d dihedral angles, we need the formula (|A.15P giving the 
contraction of two structure constants. Applying this to (15. 2p . we get 



Aij± Eik± ■ Mij-i-Eji± + Aijk±Ajii± cos pikji± + cos (j)ijk± cos (j)jii± 
cos&ikji± = , — , = ^— ^—7 7 (5.5) 



where we introduced 

Eik± ■ Mij±Eji± 
cos pik,ji± = , . . • (5.6) 

\/^ik±^jl± 

If the three triangles {ij},{ik} and {jl} share an edge we obtain from the edge simplicity constraint 
that cos Pik,ji+ = cos Pikji- and hence we have that the 4d dihedral angles are equal up to a sign and 
cos 9ikji+ = cos Oikji- on the geometric (and topological) sector. 



This is in close analogy with the state of affairs in 3d space-time, where one defines the normal to the triangle to be 
n" = e°''"'e\e'j and e" are the triad vectors. 
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6 The phase space for a 4— simplex 

Unfortunately, the system of simplicity and Gaui3 constraints is, in general, quite redundant. This 
complicates the determination of the Dirac brackets. To obtain an irreducible set of constraints, we 
will therefore consider a very simple triangulation, namely the boundary of a 4-simplex, consisting 
of five tetrahedra, where every tetrahedron is glued to the other four tetrahedra. The first part of 
the reduction process, namely the one which says that the left-handed geometry coincides with the 
right-handed geometry, carries over to general triangulations. But we will see, that we have left-over 
constraints, which we shall implement later on. 

The (3d) boundary triangulation of a 4-simplex has 10 triangles {ij},i < j and 5 tetrahedra 
{i},i = 1, ... ,5. We have therefore 6 x 10 configuration variables and 6 x 10 momentum variables. 
Since we have 6 constraints per tetrahedron due to the closure constraints, we obtain 30 (first class) 
constraints. (Although one might first think that there is some redundancy in the closure constraints, 
this is not the case. Redundancy occurs only in the case that the momentum variables M^^ are 
trivial.) Therefore, we expect 60 phase space variables in the Gaufi-reduced sector. 

These 60 variables can be chosen from the following set of 140 variables 






Aij± :-- 


- Eij± • Eij± 


Aijk± '■'- 


= Eij± ■ Eikj- 


cos%i± := 


= cos Oiiji± = 



y/Niji± ■ Niji± Njik± ■ Njii± 

The final set of variables may be somewhat puzzling at first sight, but we should notice that they are 
an artefact of the connectivity of a 4-simplex. Each edge is contained in exactly three triangles and 
three tetrahedra. 

The 20 variables Aij± do not depend on whether they are computed in the tetrahedron {i} or 
the tetrahedron {j}, that is Aij± = Aji±. A priori, we have 60 variables Aijk±, namely 2x6 per 
tetrahedron. We can use the GauB constraints for a tetrahedron to reduce the number of independent 
variables to the 2x2 non-opposite angles (/'jjfc±|j Here, we have to choose some prescription within 
the tetrahedron {i} for ji,ki with I = 1,2 , such that the resulting angles (pijjkj are non-opposite. 
To extract the other 2x4 angles in the tetrahedron {i}, one has to multiply the GauB constraints for 
this tetrahedron with any of the i?-bivectors. In other words, one uses the 2x4 equations 

/ , Eij± ■ Eik± = 2_^ Aijk± = . (6.2) 

At this stage, we have 20 area-variables and 2 x 2 x 5 = 20 angle variables. These uniquely determine 
the left- and right-handed intrinsic geometry of the five tetrahedra. (Remember that the geometry of 
a tetrahedron is determined by six quantities, for instance the length of the edges.) Note that these 
variables do not just encode the information about the bivectors E, but also part of the information 
encoded in the holonomies M. This happens because of the condition Eij = —MijEji relating the 
bivectors in neighbouring tetrahedra by a holonomy. 

We are left to choose 20 variables for the extrinsic geometry, that is 2 variables per triangle. We 
therefore choose one of the three cos Oiji± per triangle {ij} and per it-sector. This does not necessarily 
mean that the angles 9ijk± are the same (yet) for different choices of the edge in the triangle to which 



cos cj>ijk± is given in (|5.1|l 
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they are attached, just that the three angles per sector per triangle are related to each other by 
relations involving the other variables. 

As we have seen, from the simplicity, cross simplicity and edge simplicity constraints we can 
conclude that the geometric quantities as computed from the left-handed sector should coincide with 
the ones from the right-handed sector. That is, we have the 30 constraints 



(6.3) 



s., - 


- -^13+ -^ij— 1 


ijiki = 


= ^ijiki+ - Aijiki- 


Cijl = 


= COS 6ijl+ - COS 6ijl 



with the understanding that in the last line of (j6.3p . the index / depends on the other two indices i,j. 

Note that we did not take into account all the edge simplicity constraints yet. This can be seen 
from a simple counting argument. We know that the Gaufi and simplicity constraints ensure the 
geometricity of the triangulation. Thus, the reduced phase space of this triangulation should be 
parametrised by 20 variables, for example, the edge lengths and their conjugate momenta. But after 
imposing the equality of the left-handed and right-handed sectors (j6.3p . we are still left with 30 
variables. Therefore, we must find some more constraints to further reduce this phase space. 

Already at this stage, however, the set of constraints ()6.3p is a second class system with an invertible 
Dirac matrix. The reason for this is that the matrix of Poisson brackets between these constraints has a 
triangular block structure, which allows one to compute with ease the determinant of this matrix. This 
will also hold for a general 3d triangulation. With this in mind, we shall reduce our phase space in two 
stages. Initially, we will compute the Dirac bracket of the first subset of simplicity constraints (j6.3p . 
and reduce the phase space accordingly. Only afterwards will we introduce another set of constraints, 
argue that they are necessary to capture the geometric sector and reduce once again. 

6.1 Reduction by the first subset of simplicity constraints 

Here we want to compute the Poisson brackets between the subset of simplicity constraints (j6.3p in 
order to determine the Dirac bracket. As the Dirac bracket of an arbitrary phase space function 
with the constraints (j6.3p vanishes, we expect that in the reduced phase space, the sectors are set to 
be equal and that the flow generated by quantities from the right-handed sector coincides with the 
flow generated from the corresponding left-handed quantities. That is, right-handed and left-handed 
quantities cease to commute. 

The brackets among the simplicity constraints Sij and between the simplicity and cross-simplicity 
constraints vanish, as the squared areas commute with any function of the bivectors E alone. The 
only non-vanishing bracket involving the simplicity constraint Sij is with Ciji, that is, this block in 
the matrix of Poisson brackets is also diagonal. More precisely, we have 



7^ / 1 A 

{Sij, Cijl} = 2—2 ( {^ijl+ — ^ijl-) V^ijl+ + ^ijl-) I ) 



(6.4) 



where 



^, .^^ ^ Eij± ■ {Niji± X {Mij±Njii±)) ^g ^^ 

y/Nijk± ■ Nijk± Njik± ■ Njik± 

as calculated in more detail in AppendixO Note that the first summand (Sjj7+ + S.y7_) vanishes on 
the constraint hypersurface due to equation ()D.4p . 
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Next we will consider the brackets among the cross simplicity constraints. Here the cross-simplicity 
constraints associated to different tetrahedra commute, so that the only non-vanishing brackets are 
between the pairs of constraints from the same triangle. That is the matrix-block associated to the 
cross-simplicity constraints consists of 5 antisymmetric blocks of size 2 x 2 on the diagonal. We have 

{Sijk, Siji} = 2 _ 1 ( (.^ijki+ — Vijki-) (Vijki+ + Vijki-) j , (6.6) 

where 

Vijki± := Eij± ■ {Eik± X Eii±) . (6.7) 

Furthermore, up to a sign factor Vijki± does not depend on the choice of the three disjoint triangles 
{ij}, {ik} and {i/} (because of the Gaufi constraint and antisymmetry of C ) and is proportional to 
the volume squared of the tetrahedron. Note that here the second summand {Vijki+ — Vijki-) vanishes 
on the geometric subsector, see appendix [Pl 

We will not need the other Poisson brackets between the constraints for our line of arguments. 
The structure for the Dirac matrix of Poisson brackets Djk '■= {C/,Cft-}, {C/} = {Sij, Sijjkj,Ciji} 
which we determined so far is 

/ 21 \ 

(6.8) 

where 21 is a diagonal 10 x 10 matrix giving the Poisson brackets between Sij and Ciji. Moreover, 53 is 
the matrix of Poisson brackets between the cross simplicity constraints Sijjkj and has 5 antisymmetric 
2x2 blocks on its diagonal. The inverse is of the following form 

* * -2t-i \ 

* «B-i . (6.9) 
2t-i / 

This already allows us to compute some of the Dirac brackets with respect to the first subset of 
simplicity constraints defined by 

{/,5}i = {f,9} - {f,CK}{D-')KL{CL,g} . (6.10) 

A straightforward calculation gives for the Dirac brackets between areas and dihedral angles 

{Aije, cos Oil jq'^i}i = 5(ij)(^iiji^T.ijk+ , (6.11) 

where e, e = ±. Here 5(^ij)(i'j') = 1 if the triangles (ij), {i' j') coincide (irrespective of their orientation) 
and zero otherwise. From (j6.1ip we can conclude (see appendix [E|) 

{a-ij , ^j'j'fc'li = ^{ij),{i'j') ■ (6.12) 

where Uij := 2"^/^^.]^ and 9ijk = 9ijk+. 

The Dirac brackets between two non-opposite 3d dihedral angles in a tetrahedron {i} are 

1 
{Aijke-,AijU'}i = -;^Vijki+ ■ (6.13) 

The areas Aij^ still commute with the variables Amnpe- 

The results of this section apply not only to the 4-simplex but to a general 3d triangulation, as 
the structure of the matrix of Poisson brackets between the first subset of simplicity constraints does 
not change. 
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6.2 Reduction by the gluing constraints 

We did not implement the full set of simplicity constraints yet, as part of the edge simplicity con- 
straints (j3.6p are missing. Indeed, the metricity of the triangulation implies gluing constraints, that 
is constraints involving neighbouring tetrahedra, which we will discuss in this section. 

From the last section, we conclude that we can reduce our considerations to the right handed 
sector, for instance. In the following, we will therefore drop the it-index. For the example of a 4- 
simplex we are left with 10 area variables Aij, 20 3d dihedral angles 4'iji,ki and 10 4d dihedral angles 
9iji. As is pointed out in [1], these variables will not in general be consistent with the geometry of a 
triangulation. For example, consider two neighbouring tetrahedra {i} and {j}. The four areas and 
two angles per tetrahedron allow one to compute the length of the edges (for an explicit procedure see 
m). But there is no guarantee that the lengths of the three shared edges computed with respect to 
tetrahedra {i} and {j} will coincide. Indeed, we have to impose the equality of these length variables 
by further constraints. Instead of using length variables, one can choose, as in [l], the 2d angles Uijki in 
the triangle {ij} between the edges shared by the triplets of tetrahedra {i}, {j}, {k} and {i}, {j}, {I}, 
respectively. Let us denote by Uijki the angle as computed from the geometrical data associated to 
the tetrahedron {i} and by ajM the same computed from {j}. 

Then, we have the 30 gluing constraints 

cos ttijki = cos a jiki , (6.14) 

which follow from the fact that the GauB, simplicity, cross simplicity and edge simplicity constraints 
allow us to conclude that the (dual) bivectors come from the edge vectors of tetrahedra, see also the 
discussion in [1]. The aijki can be written as a function of the 3d dihedral angles (pimn 

Nijk ■ Niji cos (piki - cos (pijk cos (f)iji 

cos aijki = — , = = \ — — . (6.15) 

^Nijk ■ Nijk Niji ■ Niji sm 0^^ sm (f>iji 

From these thirty constraints only 10 are independent (if the Gaui3 constraints are implemented). 
This follows from a linearisation of the constraints (j6.14p around an equilateral configuration, for 
example. In fact, it is well known that the internal geometry of a 4-simplex is determined by its 
10 areas or 10 length variables. (This holds at least locally, as there is a discrete ambiguity in the 
transformations between areas and length variables, several length assignments can lead to equal area 
assignments, see |23]. In the following all considerations apply to local patches of phase space, where 
this transformation is unique.) 

A linearisation of these constraints (j6.14p around the equilateral configuration reveals also that it 
is possible to make a choice of one 2d angle aijki per triangle (ij) such that the resulting system is 
irreducible. The Poisson bracketqj { ,}i between the constraints (|6.14p can be computed and show 
that the system is again second class. It is, however, more enlightening to use another but (locally) 
equivalent set of constraints. Once the constraints (I6.14p are imposed, they allow one to consistently 
compute the length and therefore also the 3d dihedral angles as a function of the areaqj. Thus, we 
can replace these constraints by ten constraints of the form 

Gil = ^ijiki — fijikiiAmn) , (6.16) 



or rather the Dirac brackets with respect to the first subset of simplicity constraints (|6.3p . but we will ignore this 
terminology and in the following drop the subindex from the bracket. 

^For a non-degenerate triangulation the number of triangles should be at least equal or bigger than the number of 
edges. Every triangle has three edges, but for a proper 3d piecewise linear manifold every edge should be shared by at 
least three triangles 

14 



where the fijjkj suce functions of the squared areas. Unfortunately, we cannot give these functions 
exphcitly. The reason being, that there is not an exphcitly known expression for the length as a 
function of the areas in a 4-simplex, since in order to obtain such an expression, one has to determine 
the roots of a polynomial equation of order greater than five. We will, however, not need an explicit 
expression for the fijk in this discussion. Another disadvantage is that whereas ()6.14p is local in the 
sense that it involves only variables of two neighbouring tetrahedra we cannot make such a statement 
for the constraints (16.16^ . This does not make a big difference for the 3d boundary of a 4-simplex, 
but can lead to rather non-local constraints for more complicated triangulations. 

The advantage of the constraints (16.16P is the simple structure of their Poisson brackets. The only 
non-commuting variables appearing in (I6.16p are the Aijk, so that 

{Gii,Gjj}i = 5ij-Vii(^jj)k(^jj)jn{i.j) , (6-17) 

where 1{IJ), k{IJ),m{IJ) are determined by the choice of indices in the definition (I6.16p . Hence the 
10 X 10 Dirac matrix of Poisson brackets of constraints consists of five antisymmetric 2x2 blocks and 
is invertible. 

Since the areas commute with the constraints (j6.14p or (j6.16p . the Dirac brackets (with respect 
to the gluing constraints) involving areas do not change from the Poisson brackets { ,}i defined via 
the first subset of simplicity constraints ()6.3p . In other words, the relation (j6.1ip is ultimately still 
valid; the 4d dihedral angles are conjugate to the areas. The areas associated to different triangles still 
commute with each other. Since we can now express the 3d dihedral angles as functions of the areas 
(at least locally in phase space), these angles commute with respect to the full Dirac bracket. We are 
left with the Dirac brackets between the 4d dihedral angles. From the Jacobi identity, involving two 
dihedral angles and one area one can conclude that the Dirac bracket between two dihedral angles can 
be at most a function of the areas. In section El we will conclude by an indirect argument that the 4d 
dihedral angles are even commuting with respect to the Dirac bracket. This finishes our discussion of 
the reduced phase space associated to the boundary of a 4-simplex. 

7 The phase space for general triangulations 

Here we will consider general 3d triangulations (of the 3-sphere) and the associated reduced phase 
spaces. We will assume that these triangulations satisfy the piecewise linear manifold conditions [36] . 
As we will see most of the discussion is very similar to the one for the 4-simplex. In the end, however, 
we will find that we have additional constraints which involve the areas only and which are, moreover, 
first class. 

We will denote by Nt the number of triangles in the triangulation. Because every triangle is shared 
by two tetrahedra and every tetrahedron has four triangles, the number of tetrahedra is equal to half 
the number of triangles Nf. The number of phase space variables we start with is 2 x 6 x A'^t. We 
have 6 X -^ X Nt (first class) GauB constraints, hence the Gaufi reduced phase space has dimension 
6 X Nf. We can again choose the gauge-invariant variables from the areas, 3d dihedral angles and 4d 
dihedral angles as defined in (16. ip . More precisely we have two area variables per triangle (for the two 
±-sectors) and four independent 3d angle variables per tetrahedron, that is two per triangle. Hence 
we are again left with two 4d dihedral angles Oikji± per triangle (ij). 

Therefore as a first set of constraints, we have the subset (j6.3p . which enforce the equality of the 
quantities as computed in the different chiral sectors. The Dirac matrix of Poisson brackets between 
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the constraints has the same block structure as the one for the 4-siniplex and hence is invertible. We 
will also obtain the same Dirac brackets (j6.11|6.13p as before. 

In a second step we have to impose the gluing constraints (I6.14p . Again these ensure that the length 
can be computed consistently from the areas. Hence these constraints allow us again to express also 
the 3d dihedral angles as a function of the areas, that is to introduce the two constraints Gij, 7 = 1,2 
per tetrahedron. The Dirac matrix of these constraints has again block-diagonal form and is invertible. 
Note that the original gluing constraints (53) are local in the sense that they involve only variables of 
neighbouring tetrahedra. This is not necessarily the case for the constraints Gn. 

The Gil constraints do however not exhaust the gluing constraints (|6.14p . The reason is that for a 
general 3d triangulation there are more triangles then length variables. Hence a consistent geometry 
also implies constraints Kj{A) between the areas. The number of these area constraints is given by the 
difference of the number of triangles and the number of edges in the 3d triangulation. This suggests 
that one introduce new configuration variables: namely the constraints Kj{A) and length variables 
le{A.), where the subindex e denotes the edges of the triangulations. Note that in general there are 
different possibilities to define the length variables as functions of the areaqfl^ these choices will differ 
by terms proportional to the constraints Ki{A). One way to find a reduced phase space is to find 
functions involving also the dihedral angles 6, that are conjugated to the constraints Kj and the length 
variables le, respectively. 

To be more concrete consider the example of the boundary triangulations of two glued 4~simplices. 
This triangulation has eight tetrahedra and six vertices. To simplify notation we will switch to a vertex 



based labelling, ie. denote by f = 1,...6 the six vertices, a^i^ = \ 2-^kim-, Gkim are the area and 
dihedral angle variableqj associated to the triangle with vertices k, I and m. Assume that we started 
with two 4-simplices a^ = (1,2,3,4,5) and aQ = (1,2,3,4,6) and glued them together by identifying 
the tetrahedra rs, tq spanned by the vertices (1, 2, 3, 4) of both 4-simplices. 

After we implemented the Gaufi constraint, the first subset of simplicity constraints (16. 3p and 
the constraints Gn we are left with 16 area and 16 4d dihedral angle variables. There are only 14 
edges, hence we have two area constraints. These two left-over constraints between the areas can be 
understood in the following way. From the areas of each 4-simplex we can compute the 3d dihedral 
angles, in particular the ones in t^ and tq. Note that these dihedral angles do not count as variables 
of our phase space in the first place, as the tetrahedra T5,tq are not part of the 3d boundary of the 
triangulation. In general we will arrive at different 3d dihedral angles for these two tetrahedra. Indeed 
the constraints between the areas impose that the 3d dihedral angles in T5 and rg have to coincide. 
This will guarantee that the six edge lengths are independent of the 4-simplex (175 or (jg) in which 
they are computed. Since the geometry of a tetrahedron can be parametrised by the four areas and 
two non-opposite 3d dihedral angles, we end up with two commuting constraints between the areas 
of the form 

Ki = (l)^j{Aik5,Amno) - (pU'^ikG^Amno) , m, n, O = 1, . . . , 4 , (7.1) 

which impose that two non-opposite 3d angles (and therefore all 3d angles) in T5 and rg coincide. 

Define lm5 and lm6 , m, = 1, ... 4 to be the length between the vertices m and 5 or 6 respectively. 
This length variables can be computed unambiguously from the areas of simplex a^ or ag respectively. 
With /^„ and if^^, mn = 1, ... 4 we will denote the six lengths of the shared tetrahedron as computed 



These functions are quite complicated if one uses only areas, but it can be more easily defined using the bivectors 
directly, see appendix IFl 

^ Here we have to make one choice out of the three possibilities 15.21 for the definition of the dihedral angle. For the 
final reduced phase space this choice does not matter. 
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from the 4-simplex (T5 or 4-simplex ag. 

Since the areas a-ijk are conjugated to the dihedral angles Oij^ we will have 



Saiofi 



{lm5,Pl5} — ^ml , Pmb — 2_^ JT. 



i<j 



dai 



{i-m&,Pl&} — Sml , Pm6 — 2_^ "q7 



i<j 



dime 



Hj5 



Hj6 



5/6 T r 5/6 V^ 9aijk . , 

Pm'n'S = (>(mn),im'n') , Pmn = 2^ ^1 ^ijk , ['■^) 

i,j,k^b/b,i<j<k 

where m, n, m', n' = 1, . . . , 4. 

Now the flow of combinations of the form Pmn ~ Pmn would violate the constraints (|7.ip as it clearly 
changes the length variables as computed from the two 4-simplices in different ways. 

The 2x4 + 2x6 momenta p we defined in (j7.2p are not independent from each other. To find an 
independent set, we will construct variables conjugated to the 2 x 2 3d dihedral angles cp^ and (j)^. 

To this end consider the change of variables from the six length {Imn} of the tetrahedron rs or tq 
respectively to the four area variables and the two dihedral angles {{aijk}, {4>j }} and the associated 
Jacobi matrix. In this Jacobi matrix, we have the entries 

Hmn ■ Qij ) \' ■'^1 

where in the partial derivative the areas are kept fixed. Now we can define 



^5/6 ■^^ ■^^ uaijk dlmn 

Oimn 



m, 717^5, 6, m<n i,j,k^()/5,i<j<k 



5/6 

Then Uj will commute with the four areas of the tetrahedron t^ or tq and with the four length 

lm5 and the four length 1^6- Moreover Uj is conjugated to (j)/ . Hence a condition to exclude 
combinations of momenta whose flow would violate the two area constraints (j7.ip could be 

Fj = i/f - z/f - CI , (7.5) 

where cj is any conveniently chosen constant. Conditions such as (j7.5p . which basically take care of 
the momenta conjugated to the area constraints (j7.ip . might arise in a proper canonical analysis of 
the discretised Plebanski action (or area Regge calculus) as secondary constraints, i.e. as the result 
of taking the Poisson brackets of the area constraints with the Hamiltonian and Diffeomorphism 
constraints. As long as these conditions form a second class system with the area constraints we do 
not need their precise form as one can check that the brackets on the reduced phase space involving 
length variables and the following choice of momenta do not depend on this form. This choice of 
momenta is Pi5,Pi6,Pmn ■= Pmn + Pmn which are conjugated to Im5,lm6,^{lmn + ^mn) respectively. 
These variables commute with the area constraints and can be taken to parametrise the reduced 
phase space. 

Let us discuss shortly a general triangulation. We will have {N^ — N^,) constraints Kj{at) ,/ = 
1, . . . , (Nt — Ne) between the area variables at associated to the triangles of the triangulation. Here 
A^t, Nf, denote the number of triangles and edges respectively. On the constraint surface Kj{a) = 0, we 
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can unambiguously define length variables le as functions of the areas associated to the edges of the 
triangulation. Choose some continuation of these functions off the constraint hypersurface. Assuming 
the length variables to be independent, we can express the areas as functions of the length variables. In 
this way we can take the derivatives of the constraints Kj{at{le)) with respect to the length variables. 
Note that these derivatives vanish on the constraint hypersurface. 

Furthermore we define conjugate momenta pe associated to the edges as 

t@e 

where we sum runs over all triangles hinging on the edge e. Then the Poisson brackets of the momentum 
Pe with an area constraint Kj gives the derivative of this constraint with respect to le and hence it 
vanishes on the constraint hypersurface. Moreover pe is conjugated to le. 

A reduced phase space with respect to the area constraints can be obtained easily since both the 
length variables and the conjugate momenta commute with the area constraints (at least weakly). 
Therefore the Poisson bracket on this reduced phase space is given by the Poisson bracket between 
these quantities on the bigger phase space parametrised by areas and 4d dihedral angles. Also if one 
introduces another set of constraints (or 'gauge conditions') in order to obtain a second class set of 
constraints {Ki,Fi) the Dirac brackets between length and momenta will coincide with the Poisson 
brackets on the bigger phase space. 

In summary in the final reduced phase space we have length variables and momentum variables 
associated to the edges of the triangulation. This phase space and momenta (j7.6p can be related to a 
canonical analysis of Regge calculus [37] . 

One might wonder whether the additional conditions (j7.4p involving the 4d dihedral angles are not 
already included in the simplicity, cross simplicity and edge simplicity constraints ()3.5|3.6p we started 
with. This is actually not the case: one can construct a phase space configuration satisfying all the 
simplicity constraints with arbitrary 4d dihedral angles, which do not need to satisfy the conditions 

dzai). 

To this end consider the 3d triangulation embedded in R'^. Define edge vectors Cij = Vj — vt between 
the vertices Vi and Vj. These can be used to construct the bivectors as duals of the wedge products 
of the edge vectors. That is choose some orientation for every triangle denoted by tijk and attach 
the bivector Eijk (with the proper orientation) constructed out of the edge vectors to this oriented 
triangle. 

The holonomies Mijk are chosen to be of the form 

M,j^ = exp(A,,,C^^^e^^£;,^,) . (7.7) 

The action of these holonomies on the bivectors can be computed straightforwardly. For instance for 

(jaa'bb'cc' ^^c^^cf^^ ^^aa' cc' ^c^^c' ^ ^ ^aa'cc' ^c^^,c' ^ ^^ g^ 

with b'"" = 2e"'^'^'^e^,e'!f,b'^. Hence Mjjfc leaves the edge vectors Cjj and Cjfc invariant. In particular the 
rotation Mjjfc leaves the bivector Eijf. invariant. Therefore all the constraints (|3.5l3.6p are satisfied. 
This holds also for the Gaufi constraints as these are automatically satisfied for bivectors coming from 
a tetrahedron. Now the values for the 4d dihedral angles 9ijk will be proportional to Xijk, which can 
be chosen arbitrarily. 

18 



8 Implementing the Dynamics 

So far, we obtained only a kinematical description of simplicial geometries. One way to implement 
the dynamics is to construct Hamiltonian and diffeomorphism constraints. The crucial question is 
whether these constraints are first class. Ideally, these constraints should reproduce the space of 
solutions to Regge calculus in the following sense. Consider the 3d boundary of a 4d triangulation. 
The 4d triangulation induces canonical data on its boundary, that is, the length of the edges and the 
4d dihedral angles. If the 4d triangulation satisfies the Regge equations the associated 3d boundary 
data should satisfy the constraints. Moreover, a configuration in a gauge orbit of another configuration 
associated to a 4d solution of the Regge equations, should also be associated to a solution of the Regge 
equations. 

It is not clear if such kind of constraints exist at all. Apart from the fact that we do not have 
a closed constraint algebra yet for Regge calculus, the main problem is that a full understanding 
of diffeomorphisms for general discrete space times is missing. This also includes quantum gravity 
models involving discrete space or space times, for instance, spin foam models or canonical loop 
quantum gravity. As is pointed out in [HI [38], the understanding of diffeomorphism symmetry is 
however crucial for the discussion of anomaly freeness and the appearance of (bubble) divergences in 
these models. Hence, an understanding of these issues is crucial for the construction and evaluation 
of quantum gravity models involving discrete space times. 

Another open problem is the relation between 4d solutions of the Regge equations, the data 
these solutions induce on the boundary and the dependence on the choice of inner triangulation. In 
particular, we used so far a 'generalised boundary' [39], namely the 3-sphere, where even for classical 
and continuous general relativity the question of well-posed initial values is unclear. Furthermore, 
there is the dependence on the choice of the inner 4d triangulation: consider the 3d boundary of a 
4-simplex. A priori, this boundary could be consistent with infinitely many inner triangulations, as 
one can perform a 1 — 5 Pachner move on a 4-simplex, i.e. subdivide this 4-simplex into five 4- 
simplices. Then, the question is how many solutions to the Regge equations exist if one prescribes for 
instance the length of the edges in the 3d boundary. One solution, which we will discuss in more detail 
below is the flat one, corresponding to having only one 4-simplex. (This 4-simplex can be refined 
into more 4-simplices such that the geometry is still flat.) In fact, in this case, the flat soluiton is the 
only solution (up to a small number of non-generic solutions which appear to arise as artifacts of the 
discretisation procedure [?, benni]. But for more complicated boundaries, we cannot exclude multiple 
solutions. 

To clarify these points, it might be valuable to consider small triangulations, where one could 
obtain a definite answer. This is work in progress [37] . Here we will only note that for the simplest 
non-degenerate case of closed 3d triangulations, namely the boundary of a 4-simplex (and more 
generally, boundaries of 4d triangulations without 'inner' triangles) such constraints can be obtained, 
leading to a flat dynamics. As already noted this might not be the only possibility. The 4-simplex as 
a 4d triangulation has no inner triangles, the Regge equations of motion (being linear in the deficit 
angles attached to the inner triangles) are trivially satisfied. Every 4-simplex can be embedded into 
4d flat space time, the dihedral angles for such 'flat' 4-simplices can therefore be computed from the 
length or area variables, see for instance [19] for explicit formulae. We will denote by Qij the dihedral 
angle at the triangle {ij} as a function of the areas. Since the dynamics imposes that the 4-simplex 
is flat, we can take as constraints the 10 functions 

Dij = 9ij — @ij . (8.1) 
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(We ignore a sign ambiguity, which encodes the orientation of the 3d boundary of the 4-siniplex.) 
This strategy is very similar to ideas used in [33] for (2 + 1) dimensional gravity. The commutation 
relations between the constraints (18. ID are given by 

{A,.0„} = m„M + .(^-fii). (8.2) 

The second term in brackets vanishes because of the Schlafli identity for a 4-simplex. The Schlafli 
identity 

Y^aij5Qij = Q, (8.3) 

valid for any variations 5 on the space of flat 4-simplices, guarantees that the Regge action for a 
4-simplex X^ji^uC^ ~ ®ii) i^ ^ generating function for the dihedral angles and hence the term in 
brackets in (j8.2p vanishes. 

Also the first term on the right hand side of (j8.2p vanishes. As we did not explicitly compute 
the Dirac bracket between the dihedral angles, we will use an indirect argument. As explained in 
the Appendix |Fl one can construct combinations of the flatness constraints whose flow leaves the 
simplicity constraints invariant on the subspace of flat connections. The geometric interpretation of 
these constraints is to translate the vertices of the triangulation. Since the flatness constraints are 
first class this holds also for the vertex translation generators. Moreover, since these combinations of 
flatness constraints leave the simplicity constraints invariant, there are still first class with respect to 
the Dirac brackets in the final reduced phase space. 

In the following, we argue that the constraints (j8.ip are gauge invariant combinations of these 
flatness constraints and therefore must be first class. Then, we can conclude, that since the Dirac 
bracket between the dihedral angles Oij, 6^1 appearing in (j8.2p can be at most a function of the areas, 
it has to vanish in order to ensure that the constraints are first class. 

Firstly, the constraints ()8.ip implement also the flatness of the (geometric) 4-simplex. Secondly, 
the action of the constraints (18.11) on the areas is 



{aki,Dij} = S(^ki),(ij) ■ (8-4) 

Hence these constraints generate displacements of the vertices in a way such that only one area at a 
time changes. Constraints, which change only one length variable at a time can also be constructed 

{him, 2_^ Jyi %} ^ ^(klm){opq) ■ (8-5) 

jj '-'^opq 

The action of these constraints is to displace one of the vertices of the edge in question in the 
direction normal to the other three edges adjacent to this vertex. In this way one obtains four 
displacements per vertex corresponding to the four dimensions of the embedding space-time. This 
gives 4x5 vertex displacements for the 4-simplex. Of these 20 displacements, 10 combinations 
correspond to the 6 global rotations and 4 translations of the 4-simplex. Hence, we have the correct 
number, 10, of independent constraints. 

In summary, for the 3d boundary triangulation of a 4-simplex, we can obtain first class and even 
Abelian constraints, which generate the displacement of the vertices. 

This line of arguments can be extended to any 3d triangulation of a 3-sphere which can be taken 
as the boundary triangulation of a 4d triangulation without inner triangles. Given a 4d triangulation 
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without inner triangles, define the functions P^il) as derivatives of the Regge action (coinciding in this 
case with the Hamilton-Jacobi functional) with respect to the length variables. 

^e - 1- E «*(- - E ®tii)) = E Ir^^ - E 0^(0) > (8-6) 

'^ t a@t me ^ a@t 

where in the second sum we sum over the 4-siniplices a hinging on the triangle t, Q^ are the dihedral 
angles in a 4-siniplex a expressed as a function of the lengths. The terms with derivatives of the 
dihedral angles vanish because of the Schlafli identity. This suggests that one define momenta conjugate 
to the length variables by 

^-=E|r("-^*)' ^^-^^ 

t@e ^ 

which only differ by a term involving the areas from the momenta pe defined in ()7.6p and commute 
with the area constraints. 

The flatness constraints can be brought into a form Df. = p'^ — P^ih') prescribing the momenta p'^ 
as function of the length variables. The brackets between these constraints are given by 



{o.,o,} = te.p;.}-7(f -If) 



where the second term vanishes because the P^ come from a generating function. Also, the first term 
vanishes as the fiatness constraints are first class on the BF phase space with which we started. 

Note that such triangulations without inner triangles correspond to a 'tree diagram' truncation in 
group field theory, which is a method to include the sum over all triangulations into a path integral, 
see [38] . 

Let us count the number of variables and constraints for this type of triangulation. Note that 
these can be generated by applying 1 — 4 Pachner moves starting with the boundary triangulation of 
a 4-simplex. For every such Pachner move, we add one vertex and four edges to the triangulation. 
Hence we end up with zero physical degrees of freedom since we expect four (first class) translation- 
generating constraints per vertex, whereas every edge gives one configuration and one momentum 
variable. 

The situation changes if one considers triangulations of the 3-sphere, which cannot be understood 
as the boundary of 4d triangulations without inner triangles. Such triangulations can be obtained by 
applying 2 — 3 Pachner moves to the 3d triangulation, i.e. by subdividing two neighbouring tetrahedra 
into three tetrahedra sharing one (new) edge. In the 4d picture, this corresponds to gluing a 4-simplex 
with two tetrahedra sharing a triangle to the boundary of the triangulation. This triangle turns into 
an 'inner' triangle. 

Given such a 3d triangulation with some prescribed length variables, it cannot in general be 
embedded into flat 4d space. Having performed the 2 — 3 Pachner move, in order to keep the induced 
data 4d flat, the length of the new edge has to be equal to a certain function of the length of the 
other edges, see also [40j. Hence to ensure flatness of the geometry, one has to include a subset of 
constraints which are functions of the length only. Since a 2 — 3 Pachner move does not lead to new 
vertices we still have the same number of translation generating constraints as before. With a new 
edge introduced by a 2 — 3 Pachner move, we obtain not only a new length variable (which is flxed 
by a flatness constraint involving only the length variables) but also a new momentum variable. We 
can conclude that this momentum variable is also flxed by some flatness constraints, forming a second 
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class system with the constraint which fixes the length of the new edge. In this way we obtain again 
zero physical degrees of freedom and constraints generating translations of the vertices |lj 

Ultimately one aims of course not for a theory without degrees of freedom, but for one which 
reproduces solutions of Regge calculus, which in general involve curvature. As we have seen, to impose 
that the triangulation is 4d flat we have to introduce second class constraints fixing the lengths and 
momenta of some of the edges. Nevertheless a counting of the degrees of freedom still allows for the 
possibility to have four translation generating constraints at each vertex. The question is whether one 
can replace this second class set of constraints and the translation generating constraints by just a 
set of first class constraints associated to the vertices. This would allow for local degrees of freedom 
associated to the edges whose length was fixed in the model described before giving flat dynamics. 

Another possibility to impose a dynamics, that we will leave open for further research, is to 
introduce a discrete time evolution in the form of Pachner moves. Pachner moves on the 3d boundary 
triangulation can be understood as gluing or removing 4-simplices from the 4d triangulation, hence 
one can evolve a 3d triangulation and build up in this way a 4d spacetime. These Pachner moves 
should be implemented as transformations respecting the symplectic structures of the phase spaces the 
Pachner moves are acting on. As Pachner moves will change the number of tetrahedra, triangles and 
edges, we have rather to introduce maps between phase spaces based on different 3d triangulations. 

There is an alternative possibility, namely to only allow for combinations of Pachner moves that do 
not change the connectivity in the 3d triangulation, see [JT]. Triangulation changing Pachner moves 
and combinations which do not change the triangulations can be realised for (2+l)-dimensional space- 
times [33J. 

9 Conclusions and Discussion 

Let us summarise our line of arguments. We started with the canonical phase space for S'0(4) 
BF-theory on a simplicial (3d) triangulation and considered simplicity, cross simplicity and edge 
simplicity constraints that enforce the geometricity of the configurations. As a result, on the constraint 
hypersurface, the £'-fields are the duals of wedge products of vectors, which we consistently associated 
to the edges of the triangulation. Moreover, we proved that the constraints are sufficient to ensure 
geometricity in Section HI 

Afterwards, we reduced the BF-theory phase space to a phase space describing geometric configu- 
rations, in four steps. As a zeroth step ~ to simplify the discussion- we considered the Gaufi-reduced 
phase space. In this way, the rest of the constraints have an immediate geometric interpretation. 

In the first step we take care of the constraints (j6.3p . that guarantee that the areas, 3d dihedral 
angles and 4d dihedral angleq^ as computed from the right handed sector coincide with the ones 
computed from the left handed sector. 



^'^ Here we differ from 'topological gravity' in [171118] : to obtain zero physical degrees of freedom we have to introduce 
flatness constraints which are functions of the length only and are second class with some other part of the flatness 
constraints, hence do not generate translations of the vertices. One reason why we might differ from the conclusions in 
[171 118] is that there the Gaufi constraints for BF-theory which ensure that the tetrahedra are closed are replaced by 
Gaufi constraints for the triangles, that ensure that the triangles are closed. Yet it remains unclear what kind of gauge 
transformations these Gaufi constraints based on triangles are generating. 

^^Here we actually consider only a subset of the 4d dihedral angles: A priori we can deflne three angles per triangle 
and we have to select one of these angles for the parametrisation of the Gaufi -reduced phase space. These three angles 
per triangle coincide on the geometrical subsector. But before the implementation of the gluing constraints 16.141 the 
three angles could still be different. 



22 



In the second and third step we dealt with the gluing constraints ()6.14p . These constraints can 
be split into two parts. The first set is second class and can be used to solve for the 3d dihedral 
angles in terms of the area variables. The second set comprises of constraints between the area 
variables only and allow the consistent transformation to length variables. This second part is first 
class. Moreover, we have to get rid of the momenta conjugate to the area constraints. Therefore, we 
introduce further constraints, turning the first class set onto a second class one, such that the final 
phase space is parametrised by length variables and conjugated momenta. We conjecture that these 
further constraints might arise as secondary constraints, that is, from the Poisson brackets of the area 
constraints with the Hamiltonian and diffeomorphism constraints. 

There are three different forms of Regge calculus, the original one [22] based on length variables, 
area Regge calculus [23J, which uses areas as fundamental variables but has to impose constraints 
between the areas, and area-angle Regge calculus [Ij which is based on areas and 3d dihedral angles 
together with the gluing constraints (j6.14p . In a canonical analysis, all versions should finally lead to 
the same reduced phase space based on length variables and conjugated momenta. The kinematical 
phase spaces of these theories should correspond to the phase space after the first step for area-angle 
Regge calculus. To obtain the phase space for area Regge calculus one has to solve for the 3d dihedral 
angles and finally to impose the constraints between the areas and the conjugated conditions involving 
the 4d dihedral angles to obtain the phase space for length Regge calculus. 

Also, the first step of our reduction process leads us to a phase space which corresponds to (the 
GauB constraint reduction of) Loop Quantum Gravity restricted to an appropriate fixed triangulation. 
There the gauge group is SU{2) which corresponds either to the left handed or right handed sector. 
Performing again a Gaufi-reduction one would obtain area variables, 3d dihedral angles and 4d dihedral 
angles. Note however that we did not deal with the full set of geometricity constraints yet. Indeed if one 
considers a 'classical version' of LQG restricted to a fixed triangulation one cannot for instance assign 
length variables consistently to the edges of this triangulation - the LQG phase space is truly larger 
than one describing geometric configuration. One can understand this feature 02] as an enlargement 
of the space of configurations on which LQG quantisation is based. That is, instead of considering 
smooth Ashtekar connections (the configuration variable in LQG) one also allows for distributional 
ones. A recent discussion of this for the quantum theory can be found in [43]. Here we want to 
point out that there is an explicit set of constraints ()6.14p ensuring that length variables can be 
consistently associated to the edges. Moreover there is a subset, which is second class, hence one 
cannot impose these constraints exactly via operator equations on the LQG Hilbert space. Of course 
one can still construct semi-classical states [Hj that satisfy these constraints to the zeroth order in 
Planck's constant. 

In the Loop Quantum Gravity phase space the 3d dihedral angles, that is, quantities which one 
would associate to the intrinsic 3d geometry, do not commute with each other. However, if one 
implements the gluing constraints via Dirac's procedure one can solve for these 3d angles as functions 
of the areas and these do commute on the reduced phase space. Hence, the non-commutativity of 
intrinsic 3d quantities is related to the fact that the phase space is bigger than the one for metric 
triangulations. 

This enlargement does not appear in the continuum phase space upon which LQG is based. More 
concretely in LQG one usually starts with the cotriads Ca as one half of the phase space variables. 
These are one-forms. It is, however, not possible to obtain a quantum representation for these cotriads 
- instead one quantises the two-forms E'^j, := eij^eie^ (integrated over surfaces). In the continuum, 
this change of variables can be done without difficulties, i.e. the number of variables does not change. 
If we consider a fixed triangulation and encode the cotriads into vectors associated to the edges, a 
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change to two-forms integrated over triangles also involves a change of the number of variables and 
leads therefore to an over-parametrisation of the original phase space based on cotriads. Hence, one 
can see this enlargement as an result of discretisation. Notice also that a similar issue does not appear 
in (2 + 1) dimensional gravity. 

Our analysis opens the question of how the dynamics defined, for instance, via the Hamiltonian 
constraints in loop quantum gravity |6j or the 'master constraint' in algebraic quantum gravity 05] 
interacts with the gluing constraints. A natural condition to ask for, is that the Hamiltonian should 
leave the subsector of geometric configurations invariant at least in an approximate sense. 

Simplicity constraints are central for the construction of spin foam models. They ensure that the 
dynamics imposed is not the one for (topological) BF theory, but the one of general relativity (at 
least in the gravitational subsector). Also recent work studying the propagation of perturbations in 
spin foams [16] show that a correct implementation of the simplicity constraints is crucial in order to 
obtain the equivalence with general relativity (in this case with Regge calculus). 

On the other hand, new spin foam models [27^ \28\ [29] allowed the matching of boundary data with 
the LQG Hilbert space (restricted to fixed triangulations) [17]. Since in LQG the canonical simplicity 
constraints are not completely implemented, the question arises of how this can be reconciled with 
a proper implementation of the simplicity constraints in the path integral, see also the discussion in 
[30] . The solution could be in the dynamics, i.e. in the properties of the amplitudes, despite the fact 
that the simplicity constraints seem not to hold at the kinematical level |48] . 

This brings us to the definition of the (classical) dynamics based on the reduced phase spaces 
we constructed in our work. To the largest extent we left the dynamics, leading to gravity, open for 
future research. We proposed two different strategies. One would be to find constraints generating a 
continuous time evolution in the form of constraints generating translations. The other would be to 
find Pachner moves leading to a discrete notion of time evolution. For the former we mainly discussed 
the possibility to base the translation generating constraints on the vertices of the triangulation. An 
alternative as proposed in [18jo ^^d used in LQG, is to consider the Hamiltonian as part of the 
constraints generating translations of the tetrahedra instead. We choose the former possibility as this 
seems to agree with the counting of degrees of freedom if one wants to impose a fiat dynamics. Similar 
results can be obtained in 3d [33j. 

For a very simple subclass of triangulations, in particular for the boundary triangulation for a 4- 
simplex, one can find a first class set of translation generating constraints leading to a flat dynamics. 
But already these simple examples could be used to consider mini-superspace reductions along the 
lines of [l9] or study the algebra of Hamiltonian and diffeomorphism constraints along the lines of 



Note that [l9] proposes the use of the Loop Quantum Gravity phase space based on small (at least 
at this stage) triangulations. Since this Loop Quantum Gravity phase space is bigger than the one for 
metric triangulation^lj, one can expect additional quantum effects due to these additional degrees of 
freedom. Here, the simple example of the 4-simplex could be also used to determine whether or not 
these additional degrees of freedom are propagating. 

We did not address the most interesting question here, which is whether it is possible to define 
a first class set of constraints leading to a non-trivial dynamics with local degrees of freedom. This 
might not be possible in the end. An alternative in this case is to use Pachner moves, that is to 
implement a discrete notion of time along the lines of consistent discretisation [12j. 



^^motivated by a counting of degrees of freedom with which we do not completely agree, as explained in section [8l 
'Although this is not the case for the example considered in [49], where o 
a 3-sphere. This issue will appear, however, for the boundary of a 4-siniplex 



^^ Although this is not the case for the example considered in [49], where only two tetrahedra are used to triangulate 
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Also ideas from consistent discretisation could be used to start from the discretised action for 
Plebanski theory and to perform a canonical analysis of this action, as it could clarify some points 
left open in this work. A continuum analysis has appeared recently in [32] and it would be valuable 
to know which points have to be altered in a discrete analysis. 

A Some information on 50(4) and its algebra 

The Lie group SO (4) is the group of orthogonal 4x4 matrices of unit determinant, the group of 
rotations about a fixed point in Euclidean space. Moreover, its Lie algebra 50 (4) is the algebra of 
antisymmetric 4x4 matrices of zero trace and has elements Jaa', where a, a' = 1,2,3,4. They satisfy 
the algebra 

[Jaa' , Jbb'] = C^"" ^'^ Jed , (A.l) 

r^ f^aa'bb'cc' 1 aa'rs bb'st dd'tr dd'cc' ( s:ab a'b'cc' • xa'b' abcc' xa'b ab'cc' xab' a'bcc'\ „y,A ^abcd ■ 

the totally antisymmetric Levi-Civita tensor. 

We utilise in the paper some labour saving notation. We denote the antisymmetric combination 
of Lorentz indices by capital letters from the beginning of the Latin alphabet 

T := T'-"'"' ' , where T is an arbitrary tensor. (A. 2) 

In particular, we define 

5^« := 25|;5f,| and ^^^ - g-^. (A.3) 

On top of this, we also invoke an idiosyncratic summation convention for these indices, namely 

T{^T^ = -Tfrf . (A.4) 

The upshot of these definitions is that 

^AB^BC^^AC ^^^ j.A^^AA'j,A' _ (^5) 

Also, the Hodge dual on the algebra is defined as 

*E^ = e^^E^ . (A.6) 

Furthermore, the following identities hold 

^AA'(jA'BC ^ ^BB' (jAB'C ^ ^CC fjABC (A 7) 

(jABC ^ mAA'^BB'^C'C^A'B'C'^ (^g) 

^AA'j^jA'B ^ e^B'j^AB' ^ (^g) 

(M^i)^^ = M^^ . (A.IO) 

for M G S0{4), and we define a anti-symmetric product of bi-vectors as: 

(Ti X Ts)^ := C^^^e^^rf rf . (A.ll) 

The Lie algbera splits into selfdual and anti-selfdual parts. Explicitly, we invoke projectors 

P^^':=i(5^^'±e^^'). (A.12) 
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It is straightforward to check that these projectors are orthonormal, that is P^ P^ = ^ss'Pg ■> 
furthermore P^^ + P_ = 5 . We may then proceed to generate all manner of projected quantities 

Ei = P^^'e^' , Mi = Pi^'M^' , (A.13) 

qABC ^ pAA'^A'BC ^ pBB' (jAB'C ^ pCC qABC ^ pAA' pBB' pCC" (jA'B'C /^ j^^x 

The following identity for the SO (4) structure constants 

^ABCcfB'C ^ 2^^Pi^'pBB' _ pAB'p^By (^ -^5) 

is in close analogy to the relation g^^^g" '^ ^ = (5"<^ J*^ — S""^ S"- '') for the SO(3) structure constants. 

B The discrete BF-theory 

Our continuum action is 

S[E,A]= [{E+-*E)AF . (B.l) 
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We rewrite this in the Hamiltonian formalism 



S[U, A] = JdtJ^ d'x e,,k {uf^At + n^,F/, + A^D.Uf,) , (B.2) 

where 11 := E -\ — *E, the indices i, j, k denote components of the forms on the spatial hypersurface S, 
and eijk '■= eoijfc- The components ejjfclljfc are the canonical momenta to the Ai, while the components 
Iloi are Lagrange multipliers enforcing the flatness constraint eiji^^Fj/. = 0. Also, ^0 is a further 
Lagrange multiplier imposing the GauB constraint eijkDiUjk := eijk{dilljk + [Ai,njA;]) = 0, where [,] 
denotes the commutator on the so (4) Lie algebra. As expected in a parametrisation invariant theory, 
the Hamiltonian is just the sum of constraints. The Poisson brackets between the canonical momenta 
are 

{At{x),U%{y)} = e.^kS^^'sHS - y) . (B.3) 

The discussion of the canonical phase space for discrete BF-theory follows closely |17j . In addition, 
we have to deal with the Barbero-Immirzi parameter. We discretise our 3d spatial hypersurface S 
using a simplicial lattice A. Such a simplicial complex contains tetrahedra T, triangles t, edges e 
and vertices v. We shall denote the tetrahedra by the Latin indices i,j etc. Using this, we can pick 
out the triangle t = {^j}, with respect to which the tetrahedra i and j are adjacent. On top of this, 
we may pinpoint an edge e = {ijk} which is common to the triangles {ij} and {ik}. Furthermore, 
we define the topological dual A* of the simplicial complex to be that structure which associates a 
(3 — ?7-)-dimensional object to each n~dimensional subsimplex of A. Thus, we denote the vertices of 
A* by i, j, .. and the edges of A* by pairs {ij}. The spatial connection and its canonical momentum 
possess natural discrete analogues in this setting 

A -). Aij = In Mij := In A £ so (4) , (B.4) 

Jij 

n -)■ Uij := / n e so(4) , (b.5) 

J ij 
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where Mij is the parallel transport matrix from tetrahedron i to tetrahedron j. We may then re-express 

A^A^ = ^c^^^[m-/]^^mS^. (B.6) 

The Lagrange multipliers do not possess natural discrete analogues since we are not discretising the 
full space time manifold. We write the discrete action as 



S[Uij,Mij]= j dt 



it) , (B.7) 



where the flatness and Gauss constraints are replaced by their respective discrete counterparts 

qA ^ C^BC ^MijMjk ■ ■ ■ Mnif^ , (B.8) 

Gf = ^ni. (B.9) 

k 

Thus our fundamental phase space variables are {M^^ ,Ilf-) associated to the triangles {ij). Alterna- 
tively, we may take {M[^ , Ef-) where Ef^ = 7_. (11 — -*n)^-. They satisfy the following commutation 
relations 

{Et„ Mi^] = -^ {6^^' - Y^') {6^,5,, C^'^^M,r + 6a5,, C^'^^M^^) , (B.IO) 

{Et,.E^,} = 5,,5,,-^C7^^^(5^^'-ie^^')<. (B.ll) 

The second equation may be somewhat puzzling since the momenta commute in the continuum theory, 
but this structure is natural in a theory where the momentum space is the cotangent space to a group 
manifold, in our case S0{4). In any case, such a commutation relation is necessary to satisfy the 
Jacobi identity 

{{E^,E^}, M^^} + {{E^,M^°},E^} + {{M^^, ^^}, E^} = Q, (B.12) 

where we suppressed the subindices as we consider here only variables defined on one and the same 
triangle. Furthermore, they are subject to the conditions 

Ef, = -MJf^Ef^, (B.13) 

MfM^,"" = 5^^, (B.14) 

MfM^^"" = 5^^, (B.15) 

which encapsulate the properties that the E field when viewed from one tetrahedron is rotated when 
viewed from the other, that the parallel transport matrix M is mapped to its inverse under a change 
of orientation, and that M is orthogonal. 

C The 4d dihedral angles 

We defined the (moduli of the) 4d dihedral angles 9± in (j5.2p by 

^/^ijk± ■ ^ijk± Nja± ■ Njii± 
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We can obtain an expression for sin 9± by using (sin 0±)'^ = 1 — (cos 0±)'^. By virtue of the contraction 
of two structure constants (JA.ISP we can write 

{Nijk± X {MyNjii)±) ■ {Nijk± X {M.jNjii)^) ^ ^ {Nj^k^ ■ {MijNju)±f ,^^. 

2^ Nijkj- ■ Nijk± Njii± ■ Njii± Nijk± ■ Nijk± Njii± ■ Njii± 

where we used the notation (Ti x T2)^ = C^^^e^^T^^T^ . The left hand side of (IU:2|) is therefore 
equal to (sin^-t)^. The expression on the left hand side simplifies because of 

(%fc± X {M,jN^a)±) = 2^ {E,j± ■ {M^jNju)± Eik± - Eik± ■ {MijN^ii)± Eij±) , 

= —2 Eik± ■ {MijNjii)± Eij± , (C.3) 

since the first term in the first line of ()C.3p vanishes due to the contraction of Eij± with Njii±. Hence 

we have 

g-^2g. ^ 2^ Eij± • Eij± {Eik± ■ {MijNjii)±f 

Nijk± ■ Nijk± Njii± ■ Njii± 
An expression that will appear later on in the Poisson bracket calculations is 

Eij± ■ {Nijk± X {MijNjii)±) = C± C± Eijj.Eij^Eif.^{MijNjii)j. , 

= —2 Eij± ■ Eij± Eik± ■ {MijNjii)± , 
= -2^/^sme,k,ji± \Eij±\ \Nijk±\ \Njii±\ , (C.5) 

where \T\ = \/T • T. Note that in the last line of (IC.SP we also defined the signs for 9ij^]^i±. 

D Alternative geometricity constraints 

The simplicity constraints can be partially replaced by alternative constraints, which also distinguish 
between the topological and gravitational sector. Consider the totally antisymmetric expressions 

V{E, E', E") := C^BD^DC^A^,B^„C ^^^ ^^^^ ^,^ ^„^ _ ^ABC^A^,B^„C _ ^^^^ 

in the three bivectors E,E',E". Under the duality map E i— )■ *E the expression V is mapped to W 
and vice versa. It is straightforward to check that V vanishes for a triple of bivectors of the form 

E = *{eAe'), E' = *{eAe"), E" = *{eAe"'), (D.2) 

whereas W vanishes for 

E = *{eAe'), E' = *{e'Ae"), E" = *{e"Ae). (D.3) 

Hence we have for instance y(-Ejj+,£'jfc+,-Ej;_|_) = V{Eij-,Eik-,Eii_) on the gravitational sector. 
Note that due to the Gaufi constraint, V{Eij,Eik,Eii) does not depend (modulo a sign) on the choice 
of the three triangles of the tetrahedron. Indeed it is proportional to the square of the volume of the 
tetrahedron. 

On the gravitational sector three bivectors sharing an edge (and transported to the same coordinate 
system) are of the form (ID.2P . Hence we have the constraints 

Ef, M,fC^^^e''''EfiE';{ = E^u ■ (M^jN^u) = . (D.4) 
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if the three dual bivectors involved share an edge. Indeed these constraints are used in p3] instead of 
the edge simplicity constraints. Note that 

= Eik+ ■ {M,jNja)+ + Eik- ■ (M^jNju)- , (D.5) 

so that together with (IC.5P we have sin 6ik,ii+ = — sin 6ik,ii- on the gravitational sector (and if the 
three triangles involved in the definition of Oij^i^i share an edge). 

E Poisson brackets 

Here we will compute in more detail the Poisson brackets (j6.4|6.6p . To obtain the Poisson brackets 
between the area squared variables Aij and the cosine of the 4d dihedral angles we compute 



\/Nijk± ■ Nijk± Njii± ■ Njii± J 

2E- ■-!-■ 

'^ = {Eij± , Nijk± ■ {MijNju)±} , 

yNijk± ■ Nijk± Njii± ■ Njii± 

^ 27 Eij± ■ {N,jk± X {MijNjii)^)) 

7±1 ^JNijk± ■ N,jk± Njii± ■ Njii± ' 

, 27 



"7±1 



^ik,jl± ■ (E.l) 



Here, we used that Eij ■ Eij commutes with both Nij^ and Njn. Thus, the area squared variable Aij 
commutes with all 4d dihedral angles with the exception of the ones associated to the triangle {ij}. 
As a result, the Poisson brackets between the constraints Sij and Cikji '■= cos9ij^ki+ ~ ^^^Gik.ji- ^'I'e 



1' 



{Sij,Cik,ji] — — 4-2— — Sifcj7+ ) (E.2) 

where the equality sign ~ indicates that we evaluate the expression on the gravitational sector, where 
we have Sifc,j7+ = —'^ikji- because of equations ()C.5ID.5p . The Dirac bracket {Aij^^, cos Oik jie'}i as 
defined in the main text can be computed straightforwardly to be 

{Aije, cos 9ikjle'}l — T<ikjl+ , 

~ -2^/2 sin %,fe,+ y:4~7. (E.3) 

From this we can conclude that the Dirac bracket between aij := 2~^''^x/^ij+ and 9ikji := Oikji+ is 
given by 

{aij , eik,ji}i = 1 . (E.4) 

The Poisson brackets between two non-opposite angle variables Aijk±,Aiji± can be computed to 
{Eij± ■ Eik± , Eij± ■ Eii±} = ± Eij± ■ {Eik± X Eii±) , 

=: ± — --ryiEij±,Eik±,Eii±) . (E.5) 

7 ± 1 
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Now on the gravitational sector we have V{Eij^,Eik^, Euj^) = V{Eij^, Eik-, Eu^) =: V-ijki and there- 
fore 

27 
{Sijk-, Siji} c^i — 2^— j-^ifci • (E.6) 



This gives the following Dirac brackets 



2^ 
2 



{Aijke, Aijui]i ~ -i-Vijki ■ (E.7) 



between the variables Aij^e and Aiji based on non-opposite edges in the tetrahedron {i}. 

F Flat Dynamics in the BF theory phase space 

Adapting [1^ to our purposes we want to construct combinations of the flatness constraints (|B.8P that 
leave the simplicity constraints and therefore the sector of geometric configurations invariant. 
We will use the slightly altered flatness constraints 

pA .^ ^-AcABC^CD^^DE ^ i,/^i^)(M,,M,, • • • M„,)^^ , (F.l) 

7 

where e is the edge around which the holonomy is taken and the subindex i denotes the tetrahedron 
at which the holonomy starts and ends. For a bivector Ei^^ associated to a triangle {Im} (oriented in 
the direction of the loop) sharing the edge e the action of the constraint is given by 

{EtA''F^,]^mrK^\ (F.2) 

where the equality sign ~ indicates that we evaluate the result on the subspace of flat connections 
and 9Jt;j is (some) parallel transport from the tetrahedron {i} to the tetrahedron {/}. 

Now one can define combination of flatness constraints associated to the edges hinging at a vertex 
V, such that the corresponding action translates this vertex at least on the gravitational and flat 
subsector. 

To this end one has to find a function that gives (on the gravitational subsector) the vector e" 
associated to an edge e as a function of the bivectors from one of tetrahedrons sharing this edge. 
(Unfortunately the formulae given in [17j are incorrect. The antisymmetrisations appearing in the 
volume form proposed there lead to a vanishing expression.) We can construct this edge vector by 
calculating its direction and length. The bivectors in question are 

Sij = *(ei Aea) , ^^^ = *(e2 A 63) , En = * {e^ A d) , (F.3) 

for three faces of the tetrahedron {i}. The length of the vector ei is 



11-/ ^*^'' ■ ^*^'' 

since the numerator is (lengthx volume)'^ while the denominator is (volume)'^ . Uncovering its direction 
is a somewhat more arduous calculation. We begin by defining 

^abc ^ N^^^lN^^iNf^lx'^'e"'^'"''^' = (eieieix'^'e"'^'"''^'^ n>^< + {terms proportional to e'^n^n^} , (F.5) 
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where we have chosen an auxihary vector vector x", transverse to the tetrahedron. We then use this 
to produce 

B-'^ = (A^';l"''iV;;jV;;]f) A^^'''^' = [eieieix^'e-'^''^'''') e\e\4 + {terms proportional to e^'^e^'n'^} 

(F.6) 
Finally, the vector e\ can be picked out of S"^*^, up to a multiplicative factor, by contracting it 
appropriately with the dual bivectors *£'ij and *Eii. This results in 

a Tjbcd * rpab * rpcd 

|ei| I r j^bcd * ]7;ab * JTjcd) f ]^b' c' d' * J^qh' *j^c'd'\ 

where the denominator on the right hand side normalises the resulting vector. Note that we have its 
length and direction, we know everything we need to reconstruct the vector ei. 
With these function we define the constraints associated to a vertex 

Ha{v,i) = Y,^abcdel,F'J, (F.8) 

e&v 

where Cgj are the edge vectors expressed in the coordinates of the tetrahedron {i} and we sum over 
all edges adjacent to v and oriented towards v. The orientation of the loops on which the holonomies 
are based on and the associated edges are positive. Note that the Poisson algebra of these constraints 
closes on the hypersurface defined by the vanishing of the flatness constraints Gg' = for all edges e'. 
The constraint (|F.8|) has the following action on a bivector *-E'^f = eg^eg,^., with e, e' edges adjacent 
and pointing to v. 

{Elf,X^Ha{v, i)} ^ mt^^'^'X^^l 6,,y - 9Jt^^«'A[''e"'gl, 6,/^' . (F.9) 

This is exactly the change the bivector Ei^-i undergoes under a translation of the vertex by a vector 
A (in the coordinate system {«}). Hence the vertex translation constraints (IF.8|1 leave the subsector 
of geometric configurations invariant. Note however that the entire construction works only on the 
subspace of fiat holonomies. 

G The LQG phase space 

The phase space corresponding to loop quantum gravity restricted to a fixed 3d triangulation coincides 
with the one for 5'C/(2)-BF theory. That is, we can apply the same notations and conventions as for 
the 50(4)-BF theory, we have only to change the S0(4)"labels A,B,... to SU{2)-lahels a,b,... 
taking values 1,2,3. The Poisson brackets between the basic fields are [50j 

{E^,E^,} = le^^'Et^, (G.l) 

and so on, with e the totally antisymmetric Levi-Civita tensor. 
The GauB constraints 

E^^fc = 0' (G.2) 
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generates SU{2) gauge transformations. Gauge invariant quantities can be defined analogously to the 
SO (4) case 

A rpd Tpa 

A rpa rpa 

■^ijk ■— ^ij-'^ijk 1 

COS0..,, := ^MiB^, (G.3) 

where A''".^ := e"'^'^ E\-E'^^. A counting of degrees of freedom in the gauge invariant phase space gives 
the same result as the one for say the right handed sector of SF-theory, that is we have one area 
variable Aij per triangle, two 3d dihedral angle variables Aijk per tetrahedron and we have to choose 
one Oik^ji per triangle. Again this does not necessarily mean that the three angles associated to one 
triangle coincide, just that there are relations between these 4d angles which may involve the 3d angles. 
Also the symplectic structure will be analogous to (|6.12l6.13p . That is the areas are conjugated 
to the dihedral angles 0, which encode the extrinsic curvature. Note that here {aij,6ijk}-Lq,G = 7 as 
compared to (|6.12p . which gives {aij,9ijk} = 1- The reason is that in LQG one uses the Ashtekar 
connection in which the extrinsic curvature appears multiplicated with the Immirzi parameter 7. One 
therefore needs to rescale the angles 9 by I/7 to obtain the true values for the extrinsic curvature. 
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